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In this paper we study the entropy perturbations in N-flation by using the SJ\f formalism. We 
calculate the entropy corrections to the power spectrum of the overall curvature perturbation P(. 
We obtain an analytic form of the transfer coefficient which describes the correlation between 

' the curvature and entropy perturbations, and investigate its behavior numerically. It turns out 

I that the entropy perturbations cannot be neglected in N-flation, because the amplitude of entropy 

. components is approximately in the same order as the adiabatic one at the end of inflation T^^ ~ 

^SJ ' C'(l)- The spectral index ns is calculated and it becomes smaller after the entropy modes are taken 

into account, i.e., the spectrum becomes redder, compared to the pure adiabatic case. Finally we 
study the modifled consistency relation of N-flation, and find that the tensor-to-scalar ratio (r ~ 
) ' 0.006) is greatly suppressed by the entropy modes, compared to the pure adiabatic one (r ~ 0.017) 

' at the end of inflation. 
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I. INTRODUCTION 



Inflation is now a standard paradigm for describing the physics of the very early universe, but the microphysics 
nature of the ficld(s) responsible for inflation remains unknovifn. In the last few decades intensive effort has been 
devoted to understanding the fundamental physics of the inflation theory. For simplicity, most of studies have been 
O ' focused on the effective single scalar field model, however, in the low-energy limit of string theory, more than one 
^ , scalar fields are present and they may work cooperatively to drive the infiation, such as the assisted infiation [1]. 
^ ■ Recently, Dimopoulos et al. [2] showed that the many axion fields predicted by string vacua can be combined and 
—i[ lead to a radiatively stable inflation, called N-flation. Using the random matrix theory Easther and McAllister [3] 
showed that the mass distribution for N axion fields should be in the Marcenko-Pastur spectrum form. Further, 
' many cosmological observable imprints of N-flation have been investigated, such as the tensor-to-scalar ratio r [4], 
the non-Gaussianity parameter /nl [5, 13], and the scalar spectral index ns [6] for the pure adiabatic perturbation. 
The results show that for r and /nl the deviations from the single-field models are negligible, however, the spectral 
, index ns is smaller than the case of the single-field models. The preheating process after N-flation is numerically 
' investigated in [7] and the results show that the parametric resonance is suppressed which differs significantly from 
the single-field case. 

■ Compared with the single field model, the presence of multiple fields during inflation can lead to quite different 
^\ I inflationary dynamics [8] (see [9] for a review). In particular, multiple fields can lead to the generation of entropy 
(non-adiabatic) perturbations during infiation, which can alter the evolution of the overall curvature perturbation 
^ [10] and produce a detectable non-Gaussianity [11]. For a two-field model the entropy perturbations are investigated 
• ^ both analytically and numerically [12], however, the generalization from the two-field model to the model with a 
large number of fields is less developed. For the N-flation model the entropy perturbations have been investigated 
^ , by different approaches [13, 14]. By virtue of the 6J^ formalism [15], an analytic form of spectral index is derived 
■ - - > in [13], and similar result is obtained by the authors of [14] by using a different approach [16]. In this paper, using 
the new interpretations of the SAf formalism which are developed in [17], we calculate the entropy corrections to the 
power spectrum of the overall curvature perturbation, corresponding spectral index and the tensor-to-scalar ratio. 
Our numerical results are in agreement with the earlier results in [13, 14]. 

This paper is organized as follows. In Sec. II we briefly review the constructions of N-flation and the mass distribution 
of N axion fields, then investigate the cosmological background evolutions numerically. In Sec. Ill we study the linear 
perturbation, derive explicitly the entropy corrections to the primordial power spectrum, and then calculate the power 
spectrum, spectral index and the tensor-to-scalar ratio numerically for the N-flation model. Sec. IV is devoted to our 
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conclusions. 

II. REVIEW OF N-FLATION 

In this section we briefly review the construction of N-flation, especially focus on the quadratic potential and the 
mass spectrum. Then we investigate the cosmological background dynamics with a given mass distribution. 

A. Quadratic potential in N-flation 

Dimpopoulos et al. in [2] consider a potential of N axions as 

N N 

where Vj is the periodic potential which arises solely from non-perturbative effects, // is the axion decay constant 
and A/ is the dynamically generated scale of the axion potential that typically arises from an instanton expansion. 
This scale can be many orders of magnitude smaller than the Planck scale. 

For small field values ipi <C Mpi the periodic potentials can be Taylor expanded as 

Vi{m)-\mWi + -- - , (2) 

with the masses = Aj/ff. Consider the case in which the masses {mi} are distributed uniformly and the axion 
fields start out displaced from the minimum by y/o = ajMpi, with the maximum displacement set by each axion 
decay constant 



< > (3) 



,2 ^ // 

then it is efltoctivcly equivalent to the scenario of a single field $ with a super-Planckian displacement aMpi . This 
means that the typical initial condition in the large N limit is expected to be super-Planckian and it is suitable for 
chaotic inflation. In this sense, N-flation realizes the m?^^ inflation in a very well-controlled string theory setting. 

The authors of [2] assume a uniform axion mass spectrum for simplicity, however, for a realistic model we should 
exactly determine which sorts of mass spectra are possible in string compactification. Surprisingly, in the N ^ oo 
limit, using the random matrix theory Easther and McAllister [3] obtained an essentially universal mass spectrum, 
without invoking details of the compactification, such as the intersection numbers, the choice of fluxes, or the location 
in moduli space. 

B. Mass spectrum 
Consider the Lagrangian of axions with kinetic and potential terms as 

jC = ^M^^KuV^^'W^' -V, (4) 
where the supergravity potential and the KKLT superpotential reads [18] 

pi pi 
W = Wo{S,Xa) + J2MXa)e^p{-ai{pi-m)} , (6) 

with A, B run over the dilaton S, the complex structure moduli Xai and the Kahler moduli p/. Inserting (6) into 
(5), expanding the potential around the origin ^pi = 0, and using the F-flatness conditions DaW\^j^o = 0, we have 

V = {2nfMij^'^^ + ■■■ , (7) 
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where the mass matrix is ^ 

Mij = (k^^'DaCiDbCj - ZCiCj) . (8) 

pi 

From (4) and (8) , it is easy to see that the kinetic terms and the mass matrix are obviously not diagonal in the basis 
where the superpotential is simple. We now perform two orthogonal rotations to diagonalize Kjj and M/j. First, we 
rotate the basis to make the axion kinetic terms be canonical 

= \d^^id^^' - Mij^'ifJ , (9) 

where the mass matrix becomes 

Mij = {2nf^Of{DACMD^CL - ^CMCi^O'j . (10) 
JiJJ ^ ' 

We perform the second orthogonal rotation so that M/j is diagonalized, because the potential simply takes a purely 
quadratic form. In order to reach a more clear result, it is helpful to introduce a new {N + P) x N rectangular matrix 

Rai = 27re^/V7'0/ {DaWj) , (11) 

then the mass matrix Mjj becomes 

Mij = RiaR^ . (12) 

Because we do not know about the individual terms DaCi, in what follows we will regard them as random variables, 
i.e., we take Mu as a random matrix. Now the task is to determine the eigenvalue spectra of the N x N random 
matrix. Surprisingly, with the random matrix theory, authors of [3] find that in the large N limit the mass spectrum 
is independent of concrete values of DaCi and then obtain an essentially universal mass spectrum 

= , (13) 

for a < w? < b, where (3 = N/ {N + P) is the ratio of the dimensions of the rectangular matrix R, cr^ is the variance 
of the entries of Rai and a, b are defined as 
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a = v^) , (14) 

6 = a^(l + ^/^y . (15) 



The spectrum (13) is nothing, but the Marcenko-Pastur (MP) spectrum. The normalized MP spectrum (13) describes 
the mass distribution probability of a single axion field. On the other hand, the law of large numbers ensures that 
the mass distribution of A'' axion fields obeys the distribution probability of the single field. In practice we uniformly 
split the mass range of N axions (a, 6) into N {N <^ N) bins 

{ml, ffil), (rhf , m^), • • • , {m%_^,m%) , (16) 
where mj and the width of each bin 6 are 

ml = a, = 6 , mj = mj_i+S, 6 = {fhj^ - ml) / N . (17) 

Furthermore we set the masses of the axions in each bin as 

mf = (mf_i + mf )/2 , / = 1, 2, • • • , TV . (18) 



^ Here we emphasize that the moduh Xoi Pit which appear in Ci = Aje ^'^f, are not dynamical variables. 
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Dimensionless mass variable x 

FIG. 1: This figure illustrates the mass distribution of A'' = 1500 axions versus the dimensionless mass variable xi = mj/a, in 
the case of /3 = 1/2. The red boxes denote the quantity ^i/NS in the Jth bin, and the green smoothed curve is the normalized 
mass distribution probability 'P(xi) for a single axion field (Marcenko-Pastur distribution). Note that the law of large numbers 
ensures that the mass distribution of N axions obeys the distribution probability of a single field. 



Due to the law of large numbers we then have the foUowing relation 

il=V{mj)5^ (19) 

where #/ denotes the number of axions in the /th bin. 

Note that the constraint from the renormahzation of Newton's constant requires N P, i.e., /3 ^ 1/2 [2], in this 
paper we therefore focus on the model with (3 — 1/2. By introducing a convenient dimensionless mass parameter xj 
and corresponding mass spectrum V{x), we show in Fig. 1 the mass distribution of iV = 1500 axion fields in the case 
of /? = 1/2. The parameters xj and V{x) are defined as 

XI = e = - , (20) 

a a 



2, _ V(e-x)(x-i) 



Vix) = aV{m') = ' '\ '\ '- , 1< X < C , (21) 



with X = a^/a = 1/(1 - V/3)^ 



C. Background dynamics 

In the previous subsection we have obtained the mass spectrum for the N axion fields by virtue of the random 
matrix theory. Now we discuss the cosmological background dynamics with the MP spectrum. In a flat Friedmann- 
Robertson- Walker universe, the background dynamics is described by the set of equations 



= H\ (22) 



6 , 

i^^f = -ij, (23) 
^ I 

Cpi + m>fi + n?jLpi = , (24) 
where we have taken the units with Afpj = SttG = 1. In the slow roll region 

^ = 4-. (25) 
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one has the scaling solution as 



</3/(to) 



V.7(io)/ 



where to denotes the initial time of inflation. 



(26) 




-40 -30 -20 

Number of e-folds N 

FIG. 2: The evolutions of = 1500 axion fields versus the number of e-folds A/". From top to bottom, the evolutions of seven 
fields (fi{t; x) are plotted with different dimensionless mass x — 1.055, 1.165, 1.495, 2, 5, 10, 25, respectively. 



In the large N limit, we can deal with the masse distribution of axions by employing the MP spectrum. Consequently, 
the summation over / in the equations (22) and (23) becomes integrals over mass as 



N 



N 



H = -y / ip'^{t,m^)P{'n?)dm^ 



(27) 
(28) 



In principle one can solve the evolution equations (24), (27) and (28) for one of the N fields, such as the lightest field, 
then use the scaling solution (26) to get the solutions for other iV — 1 fields during the slow roll period. Because 
of the complication of the MP spectrum V{m^), however, it is very difficult to perform the integration analytically. 
In this paper we solve the set of background equations (22) and (24) numerically (see Fig. 2), in which we consider 
N = 1500 axion fields evolving from the equal-energy initial configurations m'j(p'j{to) — mj(/3j(io) with the vacuum 
expectation value (vev) of the lightest field, (p{tQ,mi) = 1, at the initial time to- Our results show that at the initial 
stage of inflation, only the heaviest fields (such as (p\x=25) begin rolling down the potential, after a Hubble time, 
the heaviest fields are no longer over-damped. Instead of immediately becoming under-damped and oscillating they 
remain critically damped due to the existence of the lighter fields, and the potential energy of the heavier fields is 
dissipated away before it is converted into kinetic energy. As a result, inflation is mainly sustained by the lighter 
fields at the late time and ends till the lightest field is no longer over-damped. 



III. LINEAR PERTURBATIONS 



In this section using the SJ\f formalism, we investigate the entropy perturbations during inflation at the linear 
perturbation level. In the linear cosmological perturbation theory the scalar perturbations of spacetime are usually 
parameterized as 



ds^ = -(1 + 2(j))dt'^ + 2ad.,Bdtdx' 



(1 - 2^)6^j + 2d^djE 



dx'dx^ 



(29) 



where a is the scale factor, if), B and E are four scalar perturbations. One can define two important gauge invariant 
quantities as 



C = -V-i? 



5p 



(30) 
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In the following subsections we first review the SAf formalism and then analyze the entropy perturbations during 
infiation by using the method proposed by Tye, Xu and Zhang in [17], and finally calculate the power spectrum, 
spectral index and the tensor-to-scalar ratio numerically. 



A. Brief review of the 5Af formalism 

The primordial curvature perturbation ({t, x) on large scales can be usually calculated by use of the SAf formalism 
[15]. (For a multi-field infiation model the SAf formalism is nicely reviewed in [17].) One of the essential assump- 
tions of the SAf formalism is the "separate universe assumption" [19], in which separate Hubble volume evolves like 
separate Friedmann- Robertson- Walker universe where density and pressure may take different values, but are locally 
homogeneous. Due to the different e-folding numbers between separate Hubble patches, the large scale curvature 
perturbation (^{t, x) on the uniform energy density slice can be expressed as the e-folding number difference between 
the uniform energy density slice and the unperturbed spatially flat slice at the end of inflation 

atE,^) = SM = K{viiU,^),tE)-Mf , (31) 

where A/'e('y?/(i*) x), Ie) and Aff are the number of e-folds on the uniform energy density slice and spatially flat slice 
respectively. (pi{t*,x) denotes the field configurations at the time of horizon crossing and Ie stands for the time 
at the end of inflation. In general, 6Af can be expanded, up to the second order perturbations, as 

SAf = NjQ' + ^Nij (q'QJ - (g^g^)) + • • • , (32) 

where the expansion coefficients are defined as Nj = dN/dif^ , Njj = N / dif^ d'^'' and is the perturbation of ip^ 
(30) in the spatially flat gauge. 

In the multi- field scenario, it turns out convenient to identify the effective inflaton field a as the path length of the 
trajectory in the N dimensional field space 

a{t)= / y^^ie^dt, (33) 



*• 1=1 



where the vector is defined by 



N 



a'^^^j. (34) 



Furthermore, one introduces other — 1 entropy basis vectors to form a set of orthogonal basis {e„}, where 
(n = a, s) and s denotes the N — 1 entropy fields in shorthand. Then the A'' evolution equations for the background 
fields (24) can be written as the evolution equation for the effective single field a 

a + 3H& -I- = , (35) 

where the potential gradient in the direction is 



Thus the unperturbed e-folding number in the ef direction can be expressed as 



2 



rcTE rr 

Af= I —da, (37) 



^ One can prove that the e-folding number Af in (37) is equivalent to the usual definition jV = — l/M^^ 'iZi I^* ^i/^jd.'t'i ^ ^ong as 

the potential takes the decoupled form V = Vj and all fields roll monotonically during infiation. We thank Jiajun Xu for useful 
correspondence about this point. 
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and at the linear order, 5M reads 

Ie 

where denotes the field perturbations in entropy directions e^. As pointed out in [17], the first term above comes 
from the shift in the initial value a*, and it corresponds to the adiabatic perturbation in the single field case. The 
second term arises when the uniform energy slice at the end of inflation is not orthogonal to the background trajectory. 
And the third and fourth terms are both dependent on the complete inflation trajectory after f^,, which reflects the 
fact that luidcr the entropy perturbations, the inflaton follows a new trajectory with different length (the third term) 
and also different speed (the fourth term). 

For simplicity, in this paper we ignore the contributions from the second term and rewrite the last two terms using 
some geometric tricks^, then (38) becomes 

- \ —eiQidt , (39) 
t. Jt, 

with 

K = -^ + E--'y^- (40) 

One can sec from (39) that, although there exist — 1 entropy modes, can only the one which is along the Go- 
direction seed the curvature perturbation. Therefore we can use the two-field formalism [12] to discuss the entropy 
perturbation. 



H 



5N ■ 



H 

G 



B. Observational predictions 



Now we calculate the observational predictions of N-flation, such as the scalar power spectra P,^, spectral index 
Us and the tensor-to-scalar ratio r. In order to calculate the power spectrum of the curvature perturbation, it is 
convenient to move to the momentum space. The Fourier mode of Ck(i) reads 



Using the SJ\f formalism and choosing the standard Bunch-Davies vacuum, 

(QaQs) = , 

the two-point correlation functions of the curvature perturbation can be expressed as 
(Ck.(i)Ck.(t)) = N„N,{Q-Q-') + J2NsNs'{Q'Q'') , 



{QcrQa) — ' {QsQs'} — ^ss' 2j^3 



(41) 



(42) 



{§) ^l/'^l *^^<(M^;te)e;(t,)ef(M , (43) 



where the quantities with subscript * denote the quantities take the values at horizon crossing and the slow roll 
parameters are defined by 



H 



eH 



(44) 



^ The detailed derivations can be found in the Appendix A of [17]. 
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With the help of the orthogonal relation among the entropy basis vectors 



(45) 



the second term of the right hand side of (43) can be expressed as 



(46) 



with 



(47) 



2.5 



N.=-60 
N.=-50 
N.=-40 
N.=-30 
N.=-20 




-50 -40 -30 -20 -10 
Number of e-folds N 

FIG. 3: The evolutions of the transfer coefficient T^g{t,tt) versus the number of e-folds TV. From top to bottom, the 
time evolutions of T^g(t,tt) are plotted for five different wavenumbers, which cross horizon at the number of e-folds TV, = 
—60, —50, —40, —30, —20, respectively. Our results show that the amplitudes of spectra decrease with the increasing of the 
perturbation wavenumber k, which indicates a red tilt spectrum. 

For a two-field model one can argue on a very general ground that the time dependence of curvature and entropy 
perturbations in the large-scale limit can always be described by [12] 



where 



7^ ^ aHS , S = PHS 



(T a 



(48) 



(49) 



and a, P are two time dependent dimensionless functions. Integrating (48), one can obtain a general form of the 
transfer matrix which relates the curvature and entropy perturbations generated at horizon crossing to those at 
some later time t 



1 Tns\(n 
Tss [ S 



(50) 



In this paper, TZ is denoted as ^ and the sound speed Cs = 1 because the kinetic term is canonical in the N-fiation model. 
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where 



Tss{t,U) = expy' dt'l3{t')H{t')^ , Tns{t,U) = j' dt' a{t')Tss{t' ,U)H(t') . (51) 

As we have shown in (39), for the multi- field model, can only one entropy mode, which is along the direction, 
contribute to the overall curvature perturbation if we neglect the torsion in the background trajectory. That is to say, 
wc can take S as the entropy perturbation in the direction, namely, at the linear perturbation level the multi-field 
model is effectively equivalent to the two-field model. The relation (50) is still valid in the N-flation model. 
With the above results, the power spectrum of the primordial curvature perturbation can be expressed as 



where the transfer coefficient 



-^TCs(^'^*) — 2e*jVss(t, f*) 



(52) 



(53) 



measures the contribution to the overall curvature perturbation from the entropy modes. The first constant term 
in (52) comes from the pure adiabatic perturbation at the horizon crossing, and the second term, which is time 
dependent, describes the entropy contributions to the curvature perturbation. Because of the existence of the entropy 
modes, the curvature perturbation docs not conserve after the horizon crossing, and the second term docs characterize 
the time evolution of the spectrum from the horizon crossing to the end of inflation. In Fig. 3, we show the time 
evolutions of the power spectra with different wavelengthes numerically. The results show that the amplitudes of 
spectra decrease with the increasing of the perturbation wavenumber k, which implies a red tilt spectrum. 

In order to further confirm the above analysis about the power spectrum, we calculate the spectral index explicitly. 
The spectral index turns out to be 

d\ogPAtE,U,k) dlogHl dloge, '^^^^ I + 2eM,s{tE,t, 
ns-l = ^ — ; = — — ; n — r H 



= -2e--q 



d log k d log k d log k 



dlogk 



(l + 2eNss) 



(54) 



where we have fixed t = ts and the J\fss reads 



J^ss{tE,U) = —n— = 2_^—. e;(i*) 

^ cr* '-Jt 



dt. 
-2 



We plot in Fig. 4 the scalar index ns versus the logarithm of dimensionless comoving wavenumber logfc/fco- It 
can be seen from the figure that compared to the pure adiabatic case, the index becomes smaller after including the 
entropy components. This numerical results are in agreement with the analytic ones in [13, 14]. 

Finally we discuss the modified consistency relation [9, 20] in N-flation. Because the tensor perturbation is decoupled 
from the scalar one at the linear order, the gravitational wave power spectrum is frozen-in on large scales as what 
happens in the single-field model 



Pt = Pt\* = 



8Hl 
47r2 ■ 



(56) 



We define the tensor-to-scalar ratio for a given k mode which crosses horizon at e-folding number Af* as 

r(t t) = ^^(^*^ = ^* 

= e*sm'^&{t,U) , {U<t< Ie) , 

where we have introduced a dimensionless correlation angle 6 

sin 6 



(57) 
(58) 



ns 
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log(k/ko) 



FIG. 4: The scalar index ns versus the logarithm of dimensionless comoving wavenumber logfc/fco, where fco stands for the 
comoving wavenumber of the mode which crosses horizon at the number of e-folds = —60. The solid (red) curve denotes 
the contributions to ns from the pure adiabatic component, while the dashed (green) one includes both adiabatic and entropy 
components. 



We can see from (57) that, after taking into account the entropy perturbations, the tensor-to-scalar ratio (r = esin^ 8) 
is always smaller than the one for the case of the pure adiabatic perturbation (r = e,). In Fig. 5 we show the tensor- 
to-scalar ratio of the fco mode which crosses horizon at e-folding number AC = —60. The results show that the ratio 
(r ~ 0.006) is greatly suppressed by the entropy modes at the end of inflation. 
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FIG. 5: The tensor-to-scalar ratio r for the fco mode which crosses horizon at the e-folding number M* = —60. The solid (red) 
curve denotes the ratio r = e* calculated from the pure adiabatic perturbation, while the dashed (green) one r = e» sin'^ 0(t, t») 
describes the evolution of the ratio from horizon crossing to the end of inflation. 



IV. CONCLUSION 



In this paper we studied numerically the dynamics of N-flation. At the background evolution level we investigated 
the evolution of N axions with the Marcenko-Pastur mass distribution, and we found that at the initial stage of 
inflation, only the heaviest fields begin sliding down the potential, after a Hubble time the heaviest fields are no 
longer over-damped. Instead of immediately becoming under-damped and oscillating they remain critically damped 
due to the existence of the lighter fields, and all the potential energy of the heavier fields is dissipated away before it 
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is converted into the kinetic energy. As a result the inflation is mainly sustained by the lighter fields at the late time 
and ends until the lightest field were no longer over-damped. 

At the linear perturbation level, we calculated the corrections of entropy perturbations to the power spectrum of 
the overall curvature perturbation by use of the SAf formalism. We obtained an analytic form of the transfer 
coefEcient T^g, which describes the correlation between the curvature and entropy perturbations, and investigated 
its behavior numerically. Our results show that the entropy perturbations cannot be neglected in the N-flation model, 
because the amplitude of entropy components is approximately in the same order as the adiabatic one at the end of 
inflation T^g ~ 0{1). Then we calculated the spectral index ns and found that the index becomes smaller once the 
entropy modes arc included, i.e., the spectrum becomes redder than the pure adiabatic one. Finally we studied the 
modified consistency relation for the N-flation model and found that the tensor-to-scalar ratio (r ~ 0.006) is greatly 
suppressed by the entropy modes, compared to the pure adiabatic one (r ~ 0.017) at the end of inflation. 

In this paper we only considered the entropy perturbations from the third and fourth terms in the right hand 
side of (38), which depend on the whole background trajectory in field space, while ignored the corrections to the 
curvature perturbation which would arise when the uniform energy slice at the end of inflation was not orthogonal to 
the background trajectory. In addition, the additional power in the curvature perturbation, which may be produced 
by the (p)reheating process, is also out of the discussion in this paper. 
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